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Abstract 



We report an inelastic neutron scattering study of antiferromagnetic spin 
dynamics in the Haldane chain compound Ni(C2H8N2)2N02C104 (NENP) 
under external hydrostatic pressure P = 2.5 GPa. At ambient pressure, the 
magnetic excitations in NENP are dominated by a long-lived triplet mode with 
a gap which is split by orthorhombic crystalline anisotropy into a lower dou- 
blet centered at A_l ~ 1.2meV and a singlet at Ay ~ 2.5meV. With pressure 
we observe appreciable shifts in these levels, which move to A_L(2.5GPa) ~ 
1.45 meV and A||(2.5GPa) « 2.2meV. The dispersion of these modes in the 
crystalline c-direction perpendicular to the chain was measured here for the 
first time, and can be accounted for by an interchain exchange « 3 x 10~^ J 
which changes only slightly with pressure. Since the average gap value A// ~ 
1.64 meV remains almost unchanged with P, we conclude that in NENP the 
application of external pressure does not affect the intrachain coupling J ap- 
preciably, but does produce a significant decrease of the single-ion anisotropy 
constant from D/J = 0.16(2) at ambient pressure to D/J = 0.09(7) at 
P = 2.5 GPa. 



PACS numbers: 75.10Jm, 75.40Gb, 75.50Ee. 
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I. INTRODUCTION 



The Haldane ground state with short-range correlations and a gap to magnetic ex- 
citations has now been found in several quasi-one-dimensional S=l antiferromagnets 0-§l. 
The magnetic properties of these materials are well- described by the spin Hamiltonian 

^ = e|^S,-S,+i + J'|:S,-S,+5^ + D(5f)2| , (1) 

where in addition to the intrachain and interchain couplings J and J' we have included 
a uniaxial single-site anisotropy, a common feature of these systems. Because of its finite 
spin-spin correlation length, the Haldane state is stable with respect to small perturbations 
by the interchain exchange and single-ion anisotropy, and hence can be found not only in the 
pure ID Heisenberg limit of Eq. ([l|), J' = D = 0, but for a finite range of J' and D around 
this point. If, however, the interchain coupling J' exceeds some critical value Jcrit.' quantum 
disorder is no longer favored and the system enters a long-range ordered Neel ground state 
with soft Goldstone excitations. Thus, a zero-temperature phase transition is expected to 
occur in a quasi-one-dimensional S=l antiferromagnet at J' = J'crit. if the interchain coupling 
could be varied. Such quantum phase transitions have recently attracted much theoretical 
attention [0-0. Various estimates using different model approximations J^-Q as well as 
numerical studies on finite systems have been performed, giving values for the critical 
interchain coupling in the broad range 1.3 x 10~^ < J'critl J ^ 2.6 x 10~^. On the other 
hand, solid upper and lower limits on Jcrit/ J available from experiments. It is found that 
the hexagonal antiferromagnet CsNiCls with J 2.8 meV and J' / J ^ 2 x 10^^ undergoes 



3D ordering at Tn ^ 4.8 K pf, while orthorhombic Ni(C2H8N2)2N02C104 (NENP) with 



J ^ 4 meV and J' / J ^ 8 x 10 ^ demonstrates all features of the Haldane system and 



remains disordered down to millikelvin temperatures [|T6 . 

The effects of single-ion anisotropy are less drastic than those of interchain coupling, 
but are also eventually critical. As the presence of the single-ion term in Eq. (|1]) does not 
change the basic ID nature of the model, its effects may be readily explored by numerical 
simulations on finite systems of realistic size, and have been addressed in a number of studies 



10|-|T4|. These predict that strong Ising-like anisotropy [D < 0) favors Neel ordering for 



D < — 0.25J, while "easy-plane" anisotropy {D > 0) eventually drives the system into a 
so-called "planar" phase at -D ~ J where all spins have S*^ = in the ground state. 

The currently available experimental data on the effect of nonzero J' and D on the 
Haldane spin chain have been obtained from the comparison of results from different S=l 
quasi- ID antiferromagnetic compounds. This restricts the experimentally accessible region 
of the (J', D) phase diagram discussed above to just a few points. However, the exchange and 
anisotropy parameters of the spin Hamiltonian can be controlled not only by changing the 
chemical composition of the magnetic compound, but also by tuning the properties of a single 
material. The application of hydrostatic pressure can change the interatomic separations and 
local atomic environments to which the spin Hamiltonian is highly sensitive, and therefore 
provides an opportunity for the controlled and continuous tuning of the parameters in the 
Hamiltonian, albeit in a restricted range. Here we report a study of the pressure dependence 
of the magnetic properties of NENP. We have used inelastic neutron scattering to measure 
the changes in the Haldane gap modes at T ^1.8 K when an external pressure P ~ 2.5 GPa 
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is applied. We find that both the intrachain and interchain couphng are only marginally 
affected, but the single-ion anistropy is reduced substantially, and thus hydrostatic pressure 
drives this material closer to the ideal Heisenberg model and further away from any potential 
quantum phase transitions. 



II. PROPERTIES OF NENP AT AMBIENT PRESSURE 

NENP is orthorhombic with space group Pnma. It has Ni^+ chains stretching along the b 
direction with two chains per unit cell staggered in the a direction [|l^. The ambient-pressure 
(P ^ 0.1 MPa) lattice parameters measured in our experiment at T = 1.8 K after the sample 
had been pressurized are a = 15.28 A, 6 = 10.23 A, and c = 8.096 A. In describing neutron 
scattering experiments, we refer to wavevector transfer in the corresponding reciprocal lattice 
Q = hsi* + kh* + Ic* = {h,k,l). Since the Ni^+ ions are displaced by b/2 along the 
chain, it is convenient to refer to the component of wavevector transfer along the chain as 
q = Q • (b/2) = /cTT when discussing the ID behavior of NENP. At ambient pressure the 
magnetic properties of NENP are well described by the Hamiltonian of Eq. (||) with J = 
4.0(2) meV, D/ J = 0.16(2), and an interchain coupling in the a-direction J^/J = 8 x 10"'^. 
The symmetry axis of the single-site anisotropy (2;-axis) is parallel to the chain axis. The 
spin dynamics of NENP under these conditions have been carefully studied in a number of 
neutron-scattering experiments 0,0. It was shown that the principal contribution to the 
spectral density of spin fluctuations 5(Q,u;) comes from a triplet of long-lived excitations 
above a gap that follow a dispersion relation with fundamental periodicity in q of 2tt. Most 
of the spectral weight is concentrated in the close vicinity of the antiferromagnetic point 
|g — vrl < 2-7?/^ ~ 0.37r (Here ^ ~6-8 spins is the correlation length in the Haldane state). 
Over the full range 0.37r < q < tt where they have been observed 0, these excitations are 
well- described by the single mode approximation (SMA): 

. (2) 

where 1(H) jL is the ground state energy per spin. The simplest dispersion relation that has 
the correct periodicity and adequately fits all the experimental data was found to be P] 



;itu„(g) = +t;2sin2g + Asin^l (3) 
with f = 9.7 meV, and A = 34 meV^. Instead of a degenerate triplet excitation with 



a single gap A ?a 0.41J 1.64 meV at g = tt as expected in the isotropic case |113| , p!4 
the planar anisotropy in Eq. (0) splits the triplet into two branches with Ay 2.5 meV 
and Ax ~1.2 meV for fluctuations polarized parallel and perpendicular to the chain axis, 
respectively. The lower mode is further split by a small orthorhombic anisotropy £'[(S'f )^ — 
{S\f \ with E ^ O.OIJ = 0.04 meV giving A^ ^ 1.34 meV and ^ 1.16 meV This 
latter splitting is smaller than the instrumental resolution of our current experiment, and 
its effects will therefore be neglected in our subsequent discussion. The lower doublet also 
shows measurable dispersion along the a* direction, with an effective bandwidth parameter 
!S.E ^ 0.65 meV 0. This is the basis of the above estimate of the interchain exchange J^. 
The dispersion along c* has not been measured prior to this work. 
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III. EXPERIMENTAL PROCEDURE 



Our sample was a 99% deuterated NENP single crystal of mass 0.27 g, which was a part 
of the composite sample used in Ref . . The experiments were performed on the SPINS cold 
neutron triple-axis spectrometer at the National Institute of Standards and Technology. The 
beam divergences employed were 50'//ci(A~^) — 80' around the vertically focussing pyrolytic 
graphite PG(002) monochromator, and 80' — 240' around the PG(002) analyzer. A liquid 
nitrogen cooled Be filter was placed in the scattered beam path. For constant-Q scans we 
fixed the final neutron energy Ef at 5.1 meV. The resulting full width at half maximum 
(FWHM) energy resolution for incoherent elastic scattering was Ahu = 0.27(1) meV. Scans 
with constant energy transfer hu = 2 meV were performed with Ef = 4.24 meV to optimize 
the pressure cell transmission for both incident and scattered neutrons. 

To provide high pressures, we used a clamp cell for neutron scattering supplied by Oval 
Co., Ltd., which is described in detail in Ref. |TB|]. The pressure was generated at room 



temperature by applying an external load of 2x10 N on the cap of the cylindrical sample 
chamber, or "microcell." The microcell is contained inside a barrel-shaped cylinder made of 
high-density polycrystalline AI2O3, which in turn is supported by steel rings. The microcell 
accommodated the parallelepiped-shaped NENP sample (dimensions ~ 6x6x8 mm'^), as 
well as a small NaCl platelet (dimensions ^ 4 x 4 x 2mm^) used for pressure calibration. 
The microcell was filled with the pressure-transmitting fluid fiuorinert FC-75 (available 
from 3M Chemicals, Inc.). The neutron beam reaches the sample through a 10 mm- high, 
continuous central window between the steel rings after passing through the AI2O3 pressure 
cylinder and the 5 mm-thick Al outer sleeve of the cell. The loaded cell was attached to the 
cold finger of a pumped He"^ cryostat for neutron scattering and cooled down to 1.8 K. Fig. |l| 
shows the transmission through the center of the loaded cell at T = 1.8 K for neutrons in 
the energy range 2 < E < 14 meV (a reduced beam size 3 x 3 mm^ was used for this 
measurement). The complicated structure of the curve indicates that the transmission is 
limited by Bragg diffraction from the AI2O3 cylinder. In the energy range 4 < E < 14 meV, 
the measured cell transmission is well approximated by the superposition of three Gaussians, 

T{E) = 0.085 + 0.031 ■ exp |- (^%^)'| - 0.035 ■ exp |- (f^ 

+ 0.021.exp|-(^)n , (4) 



as may be seen in Fig. |^. Scattering intensities measured at P = 2.5 GPa were subsequently 
corrected using the above expression for the energy-dependent transmission of the pressure 
cell normalized to unity for the const ant- /icj = 2 meV scan. Specifically, the actual count 
rate was divided by 



T{Ef + hu)T{Ef) 



T{huj) = ^ ~ l{)JT{Ef + huj)T{EA, (5) 

/T(6.24 meV)r(4.24 meV) ^ 



where Ef is the final neutron energy, and T{E) is given by Eq. @). The actual pressure 
applied to the sample at low temperatures in our experiment was P ~ 2.5 GPa, as determined 
from the measured change Aa/a=3.2xl0~^ in the lattice constant of the NaCl sample [jl8 . 
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For the reference measurements at ambient pressure we extracted the microcell with 
the sample from the pressure cell and placed it in a standard (ILL-type) flow cryostat. Its 
transmission for 5 meV neutrons measured the same way as above was 0.60(1), which can 
be attributed to incoherent scattering from the pressure-transmitting fluid and the NENP 
sample itself. 

The sample was mounted with its (0,1,0) and (1,0,3) reciprocal lattice directions in the 
horizontal scattering plane. With this orientation, variations of the transverse (with respect 
to the spin chains) components of the wavevector transfer Q = {h, k, I) are coupled through 
h = 1/3. However, the non-zero h had negligible effect in our experiment. First, recalling 
that inelastic neutron scattering only probes spin fluctuations polarized perpendicular to Q 
T9| , we note that the (1,0,0) direction was ~ 80° out of the scattering plane. This means 



that at least 97% of the spin fluctuations polarized along a* were always probed. Second, 
since the dispersion of the excitations along a* has period 2 in /i [0, it is quite negligible for 
h < 0.2 which is the case in present study. All constant-Q scans were performed at /c = 1, 
i.e. with wavevector transfer along the chain q = it. 



IV. RESULTS AND DISCUSSION 

Our experimental results are shown in Figs. ^|-^ The scattering intensity is shown nor- 
malized to 10 minutes counting time, although due to the low transmission of the pressure 
cell the bulk of the high-pressure data was counted three times this long to obtain adequate 
statistics. In Fig. ^ we present energy scans at Q=(5/6, 1,5/2). The large component of Q 
transverse to the chains means that more than 90% of the intensity of the longitudinally 
polarized (i.e. parallel to the chain axis) fluctuations is observed. At ambient pressure 
(O.lMPa) the corresponding gap is A||= 2.49 meV, in agreement with previous studies [0,0. 
At P = 2.5 GPa this mode shifts to lower energy, and is found at A||( 2.5 GPa) = 2.2 meV. 
The numbers quoted here are determined from the flts discussed below. In contrast, the 
energy of the transversely polarized mode, which at this Q primarily contains fluctuations po- 
larized along a*, increases under pressure from Aj^ = 1.18 meV to A^(2.5GPa) = 1.43 meV. 
This increase of the energy gap in the magnetic excitations spectrum qualitatively agrees 
with the results of recent heat capacity measurements pO|] performed at pressures up to ~ 
0.5 GPa. 

Figure |^ shows scans that probe the dispersion of the transversely polarized modes with 
wavevector transfer perpendicular to the chains q± ~ 27cl. Presence of cell material in the 
beam cause much higher inelastic background and elastic scattering intensity in the high- 
pressure data. For < / < 0.5, Q makes angles between 0° and 33° with the chain direction, 
so that between 100% and 70% of the intensity of the a*-polarized transverse fluctuations 
and up to 30% of the intensity of the longitudinal fluctuations is measured. The latter is 
seen as a small rise after the main peak in the high-pressure data where the corresponding 
mode is centered at ~ 2.2 meV. As shown in Fig. ^, a small shift of Aj^ to higher energies as / 
changes from 0.5 to is clearly observed at ambient pressure. At P = 2.5 GPa, however, this 
shift is much smaller, and barely exceeds the experimental accuracy. Finally, one constant- 
huj = 2 meV scan was performed at each pressure to probe the dispersion of the transverse 
fluctuations with q along the chains. These scans are shown in Fig. ^. 

The solid curves shown in Figs. were obtained as a result of a global flt to all the data 
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at each pressure using the resolution-convoluted theoretical expression |jT9|] for the neutron 
scattering cross-section based on the single-mode approximation (Eq. (|])) for iS""(Q, u). To 
account for the dispersion perpendicular to the chains we used the dispersion relation 



/ao o.9~ a . q , ^ , „ 1 — cos (7 cos (7 I 
^c^a(Q) = JA2 +t;2sm2g + Asm2| + (AE^)2 1 ^ , (6) 

with q± = 27tI. This was derived by adding in quadrature the energy of the Haldane modes 
given by Eq. @ with the transverse dispersion obtained for magnons in spin-wave theory. 
We note that for q close to vr, Eq. (|^) was obtained by Affleck ||^ by treating the array of 
weakly interacting Haldane chains on the basis of the nonlinear sigma model. Except for the 
common parameters Ay, A_|_, v, AE± in Eq. for each scan only the individual constant 
background and gaussian intensity and width of the diffuse tail of incoherent scattering were 
varied in the global fit. Due to the limited range of q probed, the parameter A was fixed at 
the value 34 meV^ determined previously [0. The results of the fits for the parameters in the 
dispersion relation are summarized in Table 1. The values x^=3.27 and x^(2.5GPa)=0.95 
once more indicate the validity of the single mode approximation (|^). 

It can easily be shown within the framework of perturbation theory that the single- 
ion anisotropy term D J2i{S-)'^ gives a splitting of the initially isotropic Haldane triplet such 
that the average gap value remains unchanged: |A|| + |A^ = A^. This fact can also be 
established with the symmetry arguments presented in the numerical study by Meshkov \T3\, 
where to ~ 1% accuracy the splitting was found to obey 

All =Ah + ^D , A^ = Ah-Id , (7) 

so that A|| — A± = 2D. While we found that the average gap value A^ essentially did not 
change under pressure, the splitting of the gap decreased by more than 50%. This implies 
that only D is affected by hydrostatic pressure, while J remains unchanged. Using Eq. 
we arrive at the value D{2.5 GPa) = 0.09(7)J for the anisotropy constant as compared 
to 0.16 (2) J under normal conditions. This result can be given a simple interpretation: 
since the hydrostatic pressure is isotropic, it tends to pack charge distributions inside the 
crystal more symmetrically, thereby reducing the local anisotropy. From the transverse 
dispersion parameter AE± ~ 4SV2J'J, we estimate the exchange along the c* direction 
to he J'^ ^ 3 X 10~^J and J^(2.5GPa)^ 2 x 10""^ J. This is even smaller than the value 
^ 8 X 10^"^ J in the a* direction measured previously [0. Thus J'^ appears to decrease 
with pressure, but this change is at the limit of our experimental accuracy. 



V. CONCLUSIONS 



In conclusion, changes in the spin dynamics of the Haldane gap antiferromagnet 
Ni(C2H8N2)2N02C104 (NENP) under applied hydrostatic pressure of 2.5 GPa were studied 
by inelastic neutron scattering, and another point in the (D/J,JYJ) phase diagram of the 
S=l spin system in the Haldane state has been attained experimentally. The principal effect 
of pressure was found to be a decrease by a factor 1.7 of the single- ion anisotropy constant 
in the spin Hamiltonian of the Ni^^ magnetic ions. This finding explains the increase with 
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pressure of the effective spin gap in NENP observed in tlie recent heat capacity measure- 



ments |20[) cLS being the result of the increasing energy of the lower doublet component of 
the Haldane triplet. The exchange constant in the c-direction perpendicular to the chains 
was also measured for the first time. It is smaller than the exchange along a, and seems to 
decrease slightly with pressure. Thus the overall effect of applying hydrostatic pressure in 
NENP is not to move the system towards the phase transition to the Neel ordered ground 
state or the large-D planar phase but to bring it closer to the isotropic ID limit. 
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FIGURES 



FIG. 1. Energy dependence of neutron transmission for the pressure cell loaded with our 
NENP sample at T = 1.8 K. The solid line is an approximate fit given in Eq. (^). 

FIG. 2. Constant-Q scans at Q = (.833,1,2.5) picking up both transverse and longitudinal 
fluctuations at P = 2.5 GPa (top) and P = 0.1 MPa (bottom). 

FIG. 3. Scans at different q± probing the dispersion of the transverse fluctuations with wavevec- 
tor transfer perpendicular to the chains: (a)-(c) P = 2.5 GPa; (d) and (e) P = 0.1 MPa. 

FIG. 4. Dispersion curves for the transverse polarized gap mode obtained with Eq. (|^ using 
parameters from Table 1. Open symbols and error bars were obtained from independent uncon- 
strained fits of each individual peak with the resolution convoluted lineshape. 

FIG. 5. Constant-^w scans sensitive to the velocity of excitations along the chains at 
P = 2.5 GPa (top) and P = 0.1 MPa (bottom). 
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TABLES 



TABLE L Summary of the dispersion parameters obtained from the global fits to data at two 
pressures; Ah is the average gap defined with Eq. (^). 



Pressure 


Ax 


^11 




V 


AE± 


(Pa) 


(meV) 


(meV) 


(meV) 


(meV) 


(meV) 


10^ 


1.180(4) 


2.49(1) 


1.62(1) 


9.34(20) 


0.43(3) 


2.5-109 


1.426(6) 


2.20(2) 


1.68(2) 


9.47(20) 


0.32(7) 
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